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Abstract 

Wightman function, the vacuum expectation values of the field square and the energy- 
momentum tensor are investigated for a massive scalar field with general curvature coupling 
parameter inside a wedge with two coaxial cylindrical boundaries. It is assumed that the 
field obeys Dirichlct boundary condition on bounding surfaces. The application of a variant 
of the generalized Abel-Plana formula enables to extract from the expectation values the 
contribution corresponding to the geometry of a wedge with a single shell and to present the 
interference part in terms of exponentially convergent integrals. The local properties of the 
vacuum are investigated in various asymptotic regions of the parameters. The vacuum forces 
acting on the boundaries are presented as the sum of self-action and interaction terms. It 
is shown that the interaction forces between the separate parts of the boundary are always 
attractive. The generalization to the case of a scalar field with Neumann boundary condition 
is discussed. 



PACS numbers: ll.lO.Kk, 03.70. -Fk 



1 Introduction 

The nontrivial properties of the vacuum state are among the most important predictions in 
quantum field theory. These properties are manifested in the response of the vacuum to external 
influences such as external fields. A simple model of the influence is realized by imposing 
prescribed boundary conditions on the field operator. The distortion of the spectrum for the 
zero-point fluctuations of a quantum field by these conditions results in the shifts in the vacuum 
expectation values of physical observables, such as the vacuum energy density and stresses, and 
induces vacuum forces acting on constraining boundaries. This is the well known Casimir effect 
(see [U El El Hj and references therein) . The Casimir effect is common to all systems characterized 
by fluctuating quantities and has important implications on all scales, from cosmological to 
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subnuclear. In addition to its fundamental interest this effect also plays an important role in 
the fabrication and operation of nano- and micro-scale mechanical systems and has become an 
increasingly popular topic in quantum field theory. 

An interesting topic in the investigations of the Casimir effect has always been the depen- 
dence of the physical characteristics of the vacuum on the geometry of constraining boundaries. 
Analytic results can usually be found only for highly symmetric geometries including planar, 
spherically and cylindrically symmetric boundaries. Recently exact results for the Casimir force 
in geometries of a sphere and a cylinder above a plate are obtained in O |6] (see also [7]). 
Aside from their own theoretical and experimental interest, the problems with this type of 
boundaries are useful for testing the validity of various approximations used to deal with more 
complicated geometries. In the present paper we consider a less symmetric exactly solvable 
geometry of boundaries which is a combination of a wedge with coaxial cylindrical shells. The 
Casimir effect for wedge-shaped regions is well investigated in literature [H [3 [9l [lOl [TTl [12]. 
For a conformally coupled scalar and electromagnetic fields the vacuum expectation value of 
the energy-momentum tensor inside the wedge is azimuthal symmetric. In particular, the vac- 
uum energy-momentum tensor is finite everywhere apart points on the edge. This property is 
a direct consequence of the conformal invariance in the corresponding problems and does not 
take place for a non-conformally coupled scalar field. For a scalar field with an arbitrary curva- 
ture coupling parameter satisfying Dirichlet boundary condition on the wedge sides the vacuum 
energy-momentum tensor is evaluated in |13l I14j . In addition to the azimuthal dependence 
this tensor, unlike to the case of conformally coupled fields, is also non-diagonal with nonzero 
azimuthal-radial off-diagonal component. 

The investigations of quantum effects for cylindrical boundaries have received a great deal of 
attention. In addition to traditional problems of quantum electrodynamics under the presence 
of material boundaries, the Casimir effect for cylindrical geometries is also important in the fiux 
tube models of confinement |15l 116] and for determining the structure of the vacuum state in 
interacting field theories [TTj . The calculation of the vacuum energy for the electromagnetic field 
with boundary conditions defined on a cylinder turned out to be technically a more involved 
problem than the analogous one for a sphere. First the Casimir energy of an infinite perfectly 
conducting cylindrical shell has been calculated in Ref. [18] by introducing ultraviolet cutoff 
and later the corresponding result was derived by using other methods |19[ \20\ 121] . The local 
characteristics of the corresponding electromagnetic vacuum such as energy density and vacuum 
stresses are considered in [22] for the interior and exterior regions of a conducting cylindrical 
shell, and in [23] for the region between two coaxial shells (see also |24j). The electromagnetic 
vacuum forces acting on the boundaries in the geometry of two cylinders are also considered in 
Refs. |25j . In Ref. [26] scalar vacuum densities and the zero-point energy for general Robin 
boundary condition on a cylindrical surface in arbitrary number of spacetime dimensions are 
studied for massive scalar field with general curvature coupling parameter. The corresponding 
problem for the geometry of two coaxial cylindrical shells is considered in [27]. A large number of 
papers is devoted to the investigation of the various aspects of the Casimir effect for a dielectric 
cylinder (see, for instance, [11[28] and references therein). 

In the geometry of a wedge with coaxial cylindrical boundary the modes are still factorizable 
for both scalar and electromagnetic fields and the corresponding problems are exactly solvable. 
The total Casimir energy of a semi-circular infinite cylindrical shell with perfectly conducting 
walls is considered in [29] by using the zeta function technique. For a scalar field with an 
arbitrary curvature coupling parameter obeying Dirichlet boundary condition the Wightman 
function, the vacuum expectation values of the field square and the energy-momentum tensor 
in the geometry of a wedge with an arbitrary opening angle and with a cylindrical boundary 
are investigated in |13|, 114]. The corresponding Casimir densities for the electromagnetic field 
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with perfect conductor boundary conditions on bounding surfaces are considered in |3Uj . The 
closely related problem with a cylindrical shell in the geometry of a cosmic string is discussed 
in [61\ [52] for scalar and electromagnetic fields. In both scalar and electromagnetic cases the 
application of a variant of the generalized Abel-Plana formula [24] enables to extract from the 
vacuum expectation values the parts corresponding to the geometry of a wedge without the 
cylindrical shell and to present the shell induced parts in terms of rapidly converging integrals. 
This geometry is also interesting from the point of view of general analysis for surface divergences 
in the expectation values of local physical observables for boundaries with discontinuities. The 
nonsmoothness of the boundary generates additional contributions to the heat kernel coefficients 
(see, for instance, the discussion in [281 [331 EH EH] a-nd references therein). 

In this paper we investigate one-loop vacuum quantum effects for a scalar field in the ge- 
ometry of a wedge with two coaxial cylindrical shells assuming Dirichlet boundary condition 
on bounding surfaces. This geometry generalizes various special cases previously considered 
in literature for wedge-shaped and cylindrical boundaries. In addition, we also study the role 
of nonzero mass of the field quanta. The presence of boundaries eliminates the translational 
invariance and as a result the properties of the vacuum are nonuniform. The most important 
quantities characterizing the local properties of the vacuum are the expectation values of the 
field square and the energy-momentum tensor. In addition to describing the physical structure 
of the quantum field at a given point, the energy-momentum tensor acts as the source of gravity 
in the Einstein equations. It therefore plays an important role in modelling a self-consistent 
dynamics involving the gravitational field. As the first step for the investigation of vacuum 
densities we evaluate the positive frequency Wightman function. This function gives compre- 
hensive insight into vacuum fluctuations and determines the response of a particle detector of 
the Unruh-DeWitt type. Having the vacuum energy-momentum tensor we can derive the vac- 
uum forces acting on constraining boundaries evaluating the vacuum stresses at points on the 
bounding surfaces. As we will see below, in the geometry under consideration these forces are 
position dependent on the boundary and cannot be obtained by the global method using the 
total Casimir energy (on the advantages of the local method see also |36j). In the limiting case 
from the results of the present paper the local vacuum densities are obtained for the geometry of 
a rectangular waveguide (for the local analysis of quantum fields confined in rectangular cavities 
see [MlETlEH]). 

The paper is organized as follows. The next section is devoted to the evaluation of the 
Wightman function for a massive scalar field in the region bounded by two cylindrical shells 
and by the wedge walls. This function is decomposed into three parts: the first one corresponds 
to the geometry of a wedge without cylindrical shells, the second one is induced by a single 
cylindrical shell when the second shell is absent, and the third one is induced by the presence 
of the second shell. By using the formula for the Wightman function, in section [3] the vacuum 
expectation values of the field square and the energy-momentum tensor are evaluated and their 
behavior is investigated in various asymptotic regions of the parameters. In section[5|we consider 
the vacuum forces acting on bounding surfaces. For separate boundary elements these forces 
are decomposed into self-action and interaction parts. The interaction forces are investigated 
in detail and numerical examples are presented. On the example of interaction forces we also 
demonstrate the limiting transition to the geometry of a rectangular waveguide. Finally, the 
results are summarized and discussed in section [5l 

2 Wightman function 

We consider a real scalar field if inside a wedge with opening angle (po and with two coaxial 
cylindrical shells of radii a and b, a < b (see figure [Tj). For the field with curvature coupling 



3 



parameter ^ the corresponding field equation has the form 

(V' Vi + ^R + m^) ^ {x) = 0, (1) 

where R is the curvature scalar for a [D + l)-dimensional background spacetime, Vj is the 
covariant derivative operator. For special cases of minimally and conformally coupled scalars 
one has = and ^ = = {D — 1)/4Z), respectively. Here we will assume that the background 
spacetime is flat and, hence, in Eq. ([T]) we have i? = 0. As a result the eigenmodes are 
independent of the curvature coupling parameter. However, the local characteristics of the 
vacuum such as the energy density and vacuum stresses depend on this parameter. In accordance 
with the problem symmetry we will use cylindrical coordinates (r, (p^ zi, zn), N = D — 2, and 
will assume that the field obeys Dirichlet boundary conditions on bounding surfaces: 

^p\r=j = v\^=o = ^\<p=<iyo = 0) J = a, b. (2) 

These boundary conditions modify the spectrum of the zero-point fluctuations compared with 
the case of free space and change the physical properties of the vacuum. Among the most 
important characteristics of the vacuum are the expectation values of the field square and the 
energy-momentum tensor. These expectation values can be obtained from two-point functions 
in the coincidence limit. As a two-point function here we will consider the positive frequency 
Wightman function {0\ip{x)(p{x')\0) , where |0) is the amplitude for the vacuum state. This 
function also determines the response of Unruh-DeWitt type particle detectors |39j . Here we 
consider the spatial region ^ (p ^ (po, a ^ r ^ b. The formulae for the regions r ^ a and r ^ b 
are obtained in limiting cases. Note that by using the corresponding formulae we can discuss 
various combinations of boundaries in the regions ^ p ^ pQ and 4>o ^ 4> ^ 2tt. For example, 
we can consider the situation with two cylindrical shells in the first region and without shells in 
the second one. 



t> ^yy^ 1 X 




1 /' 

/ \ 


\ 

\ 

\ 
1 





Figure 1: Geometry of a wedge with two coaxial cylindrical boundaries. 

By expanding the field operator and using the standard commutation relations, the positive 
frequency Wightman function is presented as a sum over the eigenmodes: 

(0|<^(x)<^(x')|0> = ^<^„(x)v9;(x), (3) 

o. 
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where {ipa{x), ifl^^x)} is a complete orthonormal set of positive and negative frequency so- 
lutions to the field equation satisfying boundary conditions ([2]). In the region between the 
cylindrical shells, a ^ r ^ b, the eigenfunctions are specified by the set of quantum numbers 
q; = (n, 7, k), n = 1, 2, • • • , and have the form 



ifaix) = PaQgniia, ^r) sm{qn(j)) exp (ikr|| - iivt) , (4) 



where 



to = \/7^ + f^m = |kp + m^, q = 7r/(/>o, 

and r|| = (zi, zj^), k = {ki, . . . , kjy), —oo < kj < oo. In formula we have introduced the 
notation 

gqn{ia,"fr) = Yqn{ia)Jqn{ir) - Jqn{ia)Yqn{ir) , (5) 

with Jqn{z) and Yqn{z) being the Bessel and Neumann functions. The eigenfunctions ipa{x) 
defined by dH) satisfy the boundary conditions on the inner shell and on the wedge sides. The 
eigenvalues for the quantum number 7 are quantized by boundary condition ([2]) on the surface 
r = h and are solutions of the equation 

Jqn{'ia)Yqn{lh) - Yqn{ja) Jqn{jb) = 0. (6) 

In the discussion below the corresponding positive roots we will denote by 7a = CFqn^i, I = 1,2, . . ., 
assuming that they are arranged in the ascending order, aqn^i < aqn,i+i- 

The normalization coefficient Pa in O is found from the standard orthonormality condition 
for the eigenfunctions: 

d^^W J drr J d<j) (pa{x)(Pa>ix) = —5nn'Sii'5{k - k'). (7) 

By making use of the standard integral for cylinder functions (see, for instance, |40j), one finds 



with the notation 

^"'(^) = j2( m1 ^ T ^ = Va- (9) 
•Ju{z)/Ji;{r]z) - 1 

The substitution of eigenfunctions ^ into mode-sum formula ([3]) leads to the following expres- 
sion for the positive frequency Wightman function 



'2 1^ 00 00 



f ^Afi, zg gn jz, zrj a)gqn [z, zr' /a) 

n=ii=i {2^r-W^ + kla^ 
X sm{qn(l)) sin(gn(/>') exp(ikAr|| - iu) /\t)Tql^{z)\ ^, (10) 



where Ary = ry — rjj and /S.t = t — t' . As the expressions for the eigenmodes an,i are not 
explicitly known, formula (jlOp for the Wightman function is not convenient. In addition, the 
separate terms in the sum are highly oscillatory for large values of quantum numbers. For the 
further evaluation of the summation over / we apply formula [23] 



^ V-,. ^rJ.ab, ^ Kx)dx 



-- / ' 
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dx na,qnix, 7]x) [/^(xe^^/^) + hixe-""'/^)] , (11) 
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which is a direct consequence of the generahzed Abel-Plana formula (for apphcations of the 
generahzed Abel-Plana formula in investigations of the vacuum densities in the Casimir effect 
see also [51] )■ In (HH) 

and Iqn{x), Kqn{x) are the modified Bessel functions. 

As a function h{x) in summation formula (llip we choose 

,/ \ xgqn{x,xr/a)gqn{x,xr /a) . / o / 2 i ^2~^ /^o\ 

h{x) = —i exp{-^At^/x^/a'' + k^). (13) 

\Jx^ ^ kf^a^ 

The corresponding conditions for this formula to be valid are satisfied if r + r' + | At| < 26. In 
particular, this is the case in the coincidence limit t = t' for the region under consideration. As 
a result, the Wightman function is presented in the form 

00 „ 
{OMxMx')\0) = -^^-|^^sin(gn0)sin(gn</>') / rf'^ke^'^^'-ll 

^ n=l 



/O 



J|„ (x) + Ygl (x) vr 7,^ - kl 

X Gqn{ax , rx)Gqn{ax , r' x) cosh(AiY^x2 — A;2j)|', (14) 

where h{x) is defined by (|13p and we have introduced the notation 

Gqn{x, y) = Kqn{x)Iqn{y) - Iqn{x) Kqn{y) . (15) 



In the limit b ^ oo the second term in figure braces on the right of (jl4p vanishes, whereas the 
first term does not depend on b. It follows from here that the part with the first term presents 
the Wightman function for the geometry of a wedge with a single cylindrical shell of radius a. 
The corresponding problem for a massless scalar field is investigated in [14J . For points r < b the 
second term in figure braces on the right of (I14p is finite in the coincidence limit and, hence, the 
renormalization procedure for the VEVs of the field square and the energy-momentum tensor 
is reduced to the corresponding procedure for the geometry with a single shell. In addition, in 
the coincidence limit of the arguments the x-integral in (|14p is exponentially convergent in the 
upper limit. 

In formula (jl4p . the part corresponding to the geometry with a single cylindrical shell with 
radius a can be further transformed by using the identity 

9qn\X, xr / a)gqn{x , xr / a) , i \ i l i \ I Jqn{x) 

- Jqn{xr/a)Jqn{xr /a) ^ 



xHj^^Jixr/a)Hjj^Hxr'/a), (16) 

where Hqn\x), a = 1,2, are the Hankel functions. In the corresponding integral over x with the 
second term on the right of (jl6p we rotate the integration contour by the angle tt/2 for a = 1 
and by the angle —tt/2 for a = 2. Due to the well known properties of the Hankel functions, 
under the condition r + r' — \ At\ > 2a, the integrals over the arcs of the circle with large radius 
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vanish, whereas the integrals over (0, mfcm) and {0, —iakm) cancel out. Introducing the Bessel 
modified functions one obtains 



Hx)/a /-"^ J,„{xr)J,„(xr') 



oo 



By taking into account this relation, the Wightman function is presented in the form 

00 

(0|99(x)(^(3;')|0) = {ip{x)Lp{x'))^ + {ip{x)Lp{x')) ^ - ^ ^ sm{qn(j)) sin(gn0') 



n=l 



X / Ae^^^^-il r dxx^^^^j^^^^Ggn{ax,rx)Ggn{ax,r'x) 



X cosh(AtVa;^ - A;^). (18) 

In this formula, 

00 „ 

^ n=l 

X rdxx:^2!L(^^i^exp(-zAty^^Tfci:), (19) 
is the Wightman function for the wedge without cylindrical boundaries, and 
{^{x)^{x'))^ = --^5A_^sin(gn</.)sin(gn</>')y rf^ke^'^^'-ll 



n=l 

00 



is the part of the Wightman function induced by a single cylindrical shell with radius a in the 
region r > a. Hence, the last term on the right of (jlSp is induced by the presence of the second 
shell with radius h. 

An equivalent form for the Wightman function is obtained from (|18p by using the identity 



^ njVLj^gn{ax,bx)Gqn{jx,xr) Gqn {,jx, Xr ) 



j=a,b 
Kqnibx 



^V(:rr)V(xr') - i£!iMi^^„(^^)i^^„(^/), 



with the notations na = 1, n;, = — 1, and 

Iqn{x)/Iqn{y) 



(21) 



Kqn{x)Iqn{y) - Kqn{y)Iqn{x) 



This leads to the following representation for the Wightman function 

00 

(0|(/9(x)(/?(2;')|0) = ((/9(x)(/?(xO)o + - 2D-3t^d 

n=l 

X /" Ae^*^^^" r dxx^^^^f^^^Gqn{bx,xr)Gqn{bx,xr') 

J J km \l ^ 



X cosh(AtA/a;2 - ^m)- (23) 



7 



In this formula, 



^ikArii 



{ip{x)ip{x'))^^ = - ^ ^ sm{qn4>) sm.{qn(l)') j d^ke" 

^ n=l 

X r dxx f'^"//'^,^ I,n{ xr)I,n{xr ') ^^^^(^^^^2 _ fc2^) (24) 
J km hn{hx) - 

is the part induced by a single cylindrical shell of radius h in the region r < b and the last term 
on the right is induced by the presence of the second shell. Note that formulae (j20p and ()24p 
are related by the interchange a b, In ^ Kn- For a massless scalar field these formulae are 
derived in [14 



3 VEVs of the field square and the energy-momentum tensor 
3.1 Field square 

In this section we consider the VEVs for the field square and the energy-momentum tensor in the 
region between the cylindrical shells. The VEV of the field square is obtained from the Wightman 
function in the coincidence limit of the arguments. In this limit and for points away from the 
boundaries the divergences are contained in the term ((/?(x)(/?(x'))q only. The corresponding 
renormalization procedure is realized by subtracting the part for the Minkowskian spacetime 
without boundaries. By using decompositions (jlSp and (j23p for the Wightman function and 
taking the coincidence limit of the arguments, for the renormalized VEV of the field square one 
finds 

('/'^)ren = (V'^)o,ren + + (v^^)ii', (25) 

where f = a (b) for j = b (a) and the last term on the right is given by the formula 



{ip'^)jj' = -2gAD y^sin^(gn(?!>) 

n=l 

' D-3 

dxx^x^ — w?^ ^ Qj^qn{ax,bx)Gqj^{jx,rx), (26) 



m 



with the notation 



^(D+i)/2r(p _ i)/2) • ^^'^^ 

To obtain this result we have used the formula 

where B{x,y) is the Euler beta function. In formula (j26p . the term (99^)0, ren is the renormalized 
VEV for the geometry of a wedge without cylindrical shells and the term {(p^)j is induced by 
a single cylindrical shell of radius j when the second shell is absent. Hence, the last term is 
induced by the second shell of radius / . 

The formulae for single shell terms are directly obtained from (j20p and (j24p in the coincidence 
limit. By making use of formula ()28p . in the case j = a one finds 

{ip')a = -2qAnf2sm\qncl>) H dx x [x^ - m')'^ ^^Kl^{rx\ (29) 
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and the formula for (<^^)b is obtained from here by the replacements a ^ h, I K . Note that, 
as ^j^qn{x^ y) > for X < y, the both terms (v'^)j and are negative. For points away from 

the cylindrical shells the last two terms on the right of formula (|26|) are finite. Note that both 
single shell and the second shell induced parts vanish on the wedge sides (j) = 0,(po, a < r < b. 
The part diverges on the cylindrical surface r = j with the leading term 

, ^ T{{D -l)/2) 

- -(4^)(i.+i)/2|^_^-|D-i- (30) 

for points with \r/j — 1| <C | sin0|, | sin((/)o — For points near the edges (r = j,(p = 0, 4>o) the 
leading terms in the corresponding asymptotic expansions are the same as for the geometry of a 
wedge with the opening angle 0o = tt/2. The surface divergences in the VEVs of local physical 
observables are well known in quantum field theory with boundaries and are investigated for 
various types of bulk and boundary geometries (see, for example, ^ [391 SZl SHI HH HSj). 

The term in (j25p vanishes on the shell r = j and diverges on the shell r = j' . The 

corresponding surface divergences are the same as those for a single cylindrical shell of radius 
j' . It follows from here that if we present the VEV of the field square in the form 

{^%en = {ip\,ren + Yl + ^(V''), (31) 

j=a,b 

then the interference term A((/5^) is finite everywhere. Let us consider the behavior of the 
interference part in asymptotic regions of the parameters. In the limit a ^ for fixed values r 
and b, this term vanishes as a^''. In the limit 6 ^ oo and for a massless field the interference part 
tends to zero like 1/6^+^'?"^. In the same limit under the condition mb ^ 1 the interference 
part is suppressed by the factor e~'^^^ /b^^^^^^"^ . For small values of the wedge opening angle 
one has (7^1 and, hence, the order of the modified Bessel functions in the formulae for the 
VEVs is large. By using the corresponding uniform asymptotic expansions (see, for example, 
|46j ) we can see that the main contribution comes from the term with n = 1 and from the lower 
limit of the x-integral. To the leading order for the interference term we find 

, ^ 4g(^'i)/^(a/6)^^sin^(#) 
^ ^ (27r)(^+i)/2(62 _ a^){D~i)/2 ' ^''^> 

As we see, in this limit the interference part is exponentially suppressed. For points not too 
close to the cylindrical shells, similar suppression takes place for single shell induced parts. 



3.2 Vacuum energy-momentum tensor 

The VEV for the energy-momentum tensor is obtained by using the formula 

{o\Tik\o) = hm di&^ioi^ixMx'm + 

Note that in this formula we have used the form of the metric energy-momentum tensor which 
differs from the standard one by the term which vanishes for the solutions of the field equation 
(see, for instance, [IZ]). As in the case of the field square, for points away from the boundaries the 
renormalization is realized by subtracting the part corresponding to the Minkowski spacetime 
without boundaries. By using the formulae for the Wightman function and for the VEV of the 
field square, for the renormalized VEV we obtain 

{Thren = {Tho,ren + {Th, + {T^jf, (34) 



(0|^2|o). (33) 
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where j' = a (b) for j = b (a) and the non-zero components of the last term on the right are 
given by the formulae (no summation over i) 



dx x^ ix^ 



D-3 

m ) 2 U,j^qn{ax, bx) 



X {'^tqn\^Gqn{3X,rx)] - a^'^JG qn{j X , r x)] cos{2qn(t))^ 



(35) 



(f'Ar, re sin(2gn(/)) 



n=l 



^^j,qn{o,X, bx)Gqn {jx, rx) 



dx a;^(x^ — m^) 2 
2C 



rx 



Gq„(ix,rx) + (1 - AOG'q^{jx,rx 



(36) 



with G'jj{x^y) = dyGu{x,y). In formula ([35]) we have introduced notations 



«Sf [5(y)] 

ag^b(?/)] 
4? 



(4^ - 1) 



9"(y)+(i± 



9\y) 



1 - rn^r"^ /y^ 

D - 1 
2 



+ 25'(y) 



5"(y) + -ff(y)5'(y) - ^'(y) {i ± [i - 4^(1 t i)] -^j , 

(4e - 1) [5''(y) + g\y)\ - %{y)g' {y) + ^^g\y) (4^ ± 1) 

y y^ 



(37) 
(38) 
(39) 



with g{y) = Gqn{jx,y) and in (j57j) i = 0, 3, . . . , L>. In particular, for the vacuum energy 
density and stresses along directions parallel to the cylinder axis we have the relations {TQ^j-en — 
{T^)ven = . . . = (r^)ren- This property is a direct consequence of translation invariance of the 
problem along these directions. In (j34p the term {Tf')j is induced by a single cylindrical surface 
with radius j when the second shell is absent and the term {T^)jji is induced by the presence of 
the second shell. Note that the off-diagonal component {Tl)jji vanishes on the wedge sides and 
on the cylindrical shell r = j. The formulae for the components {T!^)a are obtained from (jSSp . 
([36]) by the replacements 



Q,j^qn{aX,bx) Iqn{ax)/Kqn{ax), Gqnijx,rx) Kqn{rx). 



(40) 



The formulae for {Tf')b are obtained from the corresponding expressions for {T^)a by the re- 
placements a ^ b, I K. Single shell parts in both interior and exterior regions are in- 
vestigated in |13] [T3] for a massless scalar field. These parts diverge on the shell and for 
\r/j — 1| <C I sin(/)|, | sin(</>o — the leading term in the corresponding asymptotic expansion is 
given by the formula (no summation over i) 



D{C-^D)T{iD + l)/2) 

2D^iD+l)/2 



r - ] 



0,2,...,D. 



(41) 



-D 



For the other components to the leading order one has {T^)j ~ (^2 )i ~ k ~ j| 

As in the case of the field square, the VEV of the energy-momentum tensor can be presented 
in the form 

(OII^'^IO) = {Xho + (T'hj + ^{Th, (42) 

j=a,b 

where the surface divergences are contained in the single shell parts only and the interference 
part is finite on the shells. The explicit formula for the latter is obtained by subtracting from 
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the last term on the right (j35|) and (|36p the corresponding single sheh part. It can be checked 
that the separate terms in formulae (j3ip , ()42p satisfy the standard trace relation 

T;' = ^(e-eD)V,VV + mV, (43) 

and the continuity equation VjT^ = 0. For the geometry under consideration the latter takes 
the form 

dr{rT^)+rd^Ti = 0, (44) 
dr{rTl)+rd^T^ = TI (45) 

The behavior of the VEV for the energy-momentum tensor in the asymptotic regions of the 
parameters is investigated in the way similar to that used for the field square. In the limit a ^ 
the main contribution comes from the term with n = 1 and the interference part behaves as a?'^. 
For large values of the radius of the exterior shell, b — > oo, this part vanishes as er"^"^^ l\P^~^^l'^ 
for a massive field and like 1/5^+^''^^ for a massless one. For large values of the parameter g, 
the interference term in the VEV of the energy-momentum tensor is suppressed by the factor 

In the discussion above we have considered a model where the physical interactions are re- 
placed by the imposition of boundary conditions on the field for all modes. Of course, this is an 
idealization as real physical interactions cannot constrain all the modes of a fluctuating quan- 
tum field |9l |43| [48] . In general, the physical quantities in problems with boundary conditions 
can be classified into two main groups (see also [49j ) . The first group includes quantities which 
do not contain surface divergences. For these quantities the renormalization procedure is the 
same as in quantum field theory without boundaries and they can be evaluated by boundary 
condition calculations. The contribution of the higher modes into the boundary induced effects 
in these quantities is suppressed by the parameters already present in the idealized model. Ex- 
amples of such quantities are the vacuum densities away from boundaries and the interaction 
forces between disjoint bodies. For the quantities from the second group, such as the vacuum 
densities on the boundary and the total vacuum energy, the contribution of the arbitrary higher 
modes is dominant and they contain divergences which cannot be eliminated by the standard 
renormalization procedure of quantum field theory without boundaries. Of course, the model 
where the physical interaction is replaced by the imposition of boundary conditions on the field 
for all modes is an idealization. The appearance of divergences in the process of the evaluation 
of physical quantities of the second type indicates that more realistic physical model should be 
employed for their evaluation. In literature on the Casimir effect different field-theoretical ap- 
proaches have been discussed to extract the finite parts from the diverging quantities. However, 
in the physical interpretation of these results it should be taken into account that these terms 
are only a part of the full expression of the physical quantity and the terms which are divergent 
in the idealized model can be physically essential and their evaluation needs a more realistic 
model. It seems plausible that such effects as surface roughness, or the microstructure of the 
boundary on small scales can introduce a physical cutoff needed to produce finite values for sur- 
face quantities. Another possibility, proposed in Refs. [48j, is to replace a boundary condition by 
a renormalizable coupling between the fluctuating fleld and non-dynamical smooth background 
field representing the material (for the evaluation of the vacuum energy in smooth background 
flelds see also |50]). In this model the standard renormalization procedure of quantum fleld 
theory without boundaries provides the flnite result for the quantities which are divergent in 
the boundary condition limit. An alternative mechanism for introducing a cutoff which removes 
singular behavior on boundaries is to allow the position of the boundary to undergo quantum 
fluctuations [5T|. Such fluctuations smear out the contribution of the high frequency modes 
without the need to introduce an explicit high frequency cutoff. 
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The main subject of the present paper is the investigation of the VEVs for the field square and 
the energy-momentum tensor at points away from the boundaries and the vacuum interaction 
forces between separate parts of boundaries. In the scheme where a cutoff function is used 
instead of point-sphtting, these quantities are cutoff independent and faU into the first group. 
They do not contain surface divergences and are completely determined within the framework of 
standard procedure of quantum field theory without boundaries. We expect that similar results 
would be obtained in the model where instead of externally imposed boundary condition the 
fluctuating field is coupled to a smooth background potential that implements the boundary 
condition in a certain limit [48] . 



4 Vacuum interaction forces 



In this section we investigate the vacuum forces acting on the bounding surfaces due to the 
presence of the second cylindrical shell. First of all let us consider the forces acting on the 
wedge sides. These forces are determined by the ^-component of the energy-momentum tensor 
evaluated for (j) = 0,(/>o. Note that the off-diagonal components (T'f)j and {T^)jji vanish on the 
wedge sides and, hence do not contribute to the force. The corresponding effective pressure is 
presented in the form 

P2 = P2,wedge + P2,cyl, (46) 

where P2,wedge is the vacuum effective pressure on the wedge side when the cylindrical shells are 
absent and the part P2,cyi is induced by the shells. For a conformally coupled massless scalar in 
D = 3 one has 

1 

4807rV 

The corresponding force is attractive for (pQ < n and repulsive for 0o > tt. The second term on 
the right of (|36|) is decomposed as 



P2,cyl=pi%i+pi'Sv (48) 

where P2''cyi ~ ~(^2^)jU=o is the effective azimuthal pressure on the wedges induced by a single 

cylindrical boundary with radius j, j = a, b, andp2"';?yj = — (72^)jj'|</'=o is induced by the presence 
of the second cylindrical boundary. Substituting i = 2 and (j) = 0,(pQ in the formulae for the 
VEVs of the energy-momentum tensor from the previous section, for the forces induced by the 
shells we find 

,3 



Piiy\ = - ^ / dxx{x^-w?) 2 Qj^gn{ax,hx)Ggn{jx,rx). (50) 



The expression for P^l-yi is obtained from (j49p by the replacements a ^ h, I ^ K . Single shell 
parts in the forces acting on the wedge sides, P2cyi' finite for all values r except the points 

in') 

on the edge r = j. The second shell-induced part, cyi' i^ finite for all r except the points on 

the edge r = j' , j' = a, b, j' ^ j. Note that P2"'cyi = for r = j. The integrands in (|19|) and ([50]) 
are positive and, hence, the corresponding vacuum forces are attractive. As before we can write 

?'2,cyl = P^ilyX + Ap2,cyl, (51) 

j=a,b 
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where the interference part Ap2,cyi is finite for all values o ^ r ^ 6. As it follows from 
(j50p . the corresponding forces do not depend on the curvature coupling parameter. 

In the limit a — > the main contribution into P2"cyl ^P2,cyl comes from the term with 
n = 1 and these quantities behave like a^''. In the limit 6 — > oo and for a massive scalar field the 
parts P^2lji ^^'^ ^^2,07! are exponentially suppressed by the factor e~'^^^. In the same limit and 
for a massless field the main contribution comes from the summand with n = 1 and these parts 
behave as l/b^^'^'^~^ . Now we consider the forces acting on the wedge sides in the limit of small 
values of the opening angle when the parameter q is large, (7^1. In this limit the order of the 
modified Bessel functions is large and we can use the uniform asymptotic expansions for these 
functions. By using these expansions, it can be seen that the main contribution comes from the 
n = 1 term and from the lower limit of the integral. To the leading order we find 

m g(^+^)/^exp[-2g|ln(j7r)|] 

(2vr)(^+i)/2r2|r2 - j2|(i)-i)/2- ^''^> 

In the similar way, for the interference part of the force one has: 

^P2,cyl (27r)(^+i)/2r2(52 _ a2)(^-i)/2 • ^^"^^ 

In figure [2] we have plotted the quantities a^P2"'cyi' 3 ~ ^'^'-^ o!^P2,cy\ as functions of r/a for 
D = 3 massless scalar field. The graphs are given for the wedges with 0o = 7r/2 (full curves) 
and (pQ = 37r/2 (dashed curves) and for h/a = 1.5. 




Figure 2: Vacuum forces acting on the wedge sides due to the presence of cylindrical shells for 
D = ?> massless scalar field: a^pi^ji and a^P2,cy\- The full (dashed) curves correspond to the 
wedge with (j)Q = ti/2 {(pQ = 2>it/2) and we have taken b/a = 1.5. The curves a (b) correspond 

to the effective pressure a^p^2cj\ i^'^pflyi) '^li^n the shell with radius a (6) is present only, and 
the curves ab are for the effective pressure a^P2,cyi when the both shells are present. 

Now we turn to the interaction forces acting on the cylindrical boundaries. These forces are 
determined by the ^-component of the energy-momentum tensor evaluated on the corresponding 
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surfaces. Similar to the previous case, the effective pressure on the cylindrical shell r = j is 
presented as the sum 

pii) = p^i) ^ pin') ^ (54) 

where = -{{TDq + {Tl)j)\r=j is the radial vacuum stress on the cylinder with the radius 
j when the second cylinder is absent and pi^^') = —{T^)jj/\r=j is the additional stress on this 
cylindrical surface when the second cylinder is present. Note that the off-diagonal component 
{T^)jj/ vanishes on the shell r = j and does not contribute to the force. The part p^^^ includes 
the self-action force on the cylindrical shell and belongs to the second group of quantities in the 
classification given in the previous section. Its evaluation requires more realistic model for the 
interaction of the quantum field. Unlike to the self-action force, the interaction force given by 
the second term on the right of (I54p is finite for all nonzero distances between the shells and 
can be evaluated by boundary condition calculations. From the last term on the right of (j35p 
taking i = 1 and r = j one finds: 

pin) — — siv? {qncj)) / dx x [x^ — w?') ^ ^lj^qn{ax,bx). (55) 

J „=i J"^ 

From this formula we see that pi^^') < and the corresponding forces are always attractive. 
The expression for the interaction forces between the cylindrical shells can also be written in 
the form 



pijj ) — S—i — s\v?{qn(f>) [ dxx (x^ — m^) ^ In 



Iqn{ax)Kqn{hx) 
Iqn{hx)Kqn{ax) 



(56) 



where, as before, Ua = 1, re^ = —1. As for the forces acting on the wedge sides, the interaction 
forces do not depend on the curvature coupling parameter. 

Now we consider the behavior of the interaction forces in asymptotic regions of the parame- 
ters. In the limit a — > the main contribution in the sum of formula (j55p comes from the n = 1 
term and j'^pi^^ ) ~ a^*^. For large values of the exterior shell radius, 6 — > oo, and for a massive 
field the interaction forces p'^^^ ^ are suppressed by the factor e"^™^. In the same limit and for 
a massless field one has j'^pi^i ) ~ l/h^^'^'^^^ . For small values of the wedge opening angle, 
assuming that g ^ 1, in the way similar to that used for the estimation of the forces acting on 
the wedge sides, one finds 

^ 4g(^+3)/2(,/^)2,,,,2(^^) 

In figure [3] we have plotted the interaction forces acting on cylindrical shells, a^p^^^ \ as functions 
of cj)/(l)Q for wedges with 0o = 7r/2 (full curves) and cpo = 3tt/2 (dashed curves) and for b/a = 1.5 
in the case of D = 3 massless scalar field. The curves a are for and the curves b are for 

piba) ^ 

Note that in the geometry of two coaxial cylindrical shells without a wedge the corresponding 
interaction forces are given by the formula [27j 

A ^ / f'^ D—3 

pin') = 1 ^ / duu{u^ -m^)~ flj^niau, bu), (58) 



where the prime on the sum sign means that the term n = should be halved. For D = 3 
massless scalar field and for b/a = 1.5 from this formula we have p'^'^^^ ~ — 0.437/a^ and 
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0.1 0.2 0.3 0.4 0.5 



Figure 3: Vacuum forces acting on the cylindrical shell due to presence of the second shell, a^p^^^'^ 
as functions on (j)/(j)o, for D = 3 massless scalar field. The full (dashed) curves correspond to 
the wedge with 0o = "^72 ((/"o = 37r/2) and in both cases b/a = 1.5. The curves a (b) correspond 
to the forces acting on the shell with radius a (b). 



— 0.254/a^. As it has been shown in [27], the interaction forces (j58p can also be obtained from 
the corresponding part in the total Casimir energy differentiating over the radii of cylindrical 
shells. In the geometry under consideration in the present paper the Casimir forces are position 
dependent on the boundary and cannot be obtained by global methods using the total Casimir 
energy. 

In the limit 0o ^ 0, a, 6 ^ cxd, assuming that b — a = Li and a0o = ^2 are fixed, from 
the formulae given above we obtain the corresponding results for the geometry of a rectangular 
waveguide with sides Li and L2. Here we discuss this limiting transition for the case of the 
interaction forces p^^^'\ The consideration of the other quantities is done in the similar way. In 
the limit under consideration the parameter q is large and we can replace the modified Bessel 
functions by the corresponding uniform asymptotic expansions. By using these expansions it 
can be seen that to the leading order we have 

^U^vz, « _ y = Qn. (59) 

Introducing in (j57p a new integration variable z = x/qn and by making use of ()59p . after some 
transformations, to the leading order we find 



where y = a(j). The expression on the right of this formula is the vacuum interaction force per 
unit surface between the facets of the rectangular parallelepiped separated by the distance Li 
and y is the Cartesian coordinate parallel to these facets. Other facets of the parallelepiped are 
located at y = and y = L2. Introducing in (160]) y = y' + L2/2 and taking the limit L2 — > 00 
with fixed value y', from (j60p the vacuum forces for two infinite parallel Dirichlet plates are 
obtained. Note that the local vacuum densities for a quantum field confined within rectangular 
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cavities are investigated in [SHI EH ES] (for corresponding global quantities such as the total 
Casimir energy see [H H] and references therein) . 

5 Conclusion 

In this paper we have considered one-loop quantum vacuum effects for a massive scalar field in the 
geometry of a wedge with two coaxial cylindrical shells. We have assumed that the field satisfies 
Dirichlet boundary condition on the bounding surfaces. This geometry generalizes various special 
cases previously discussed in literature, including wedge-shaped regions, cylindrical boundaries, 
and rectangular waveguides. The most important local characteristics of the quantum vacuum 
are the VEVs for the field square and the energy-momentum tensor. To evaluate these VEVs, 
as the first step we construct the positive frequency Wightman function. The corresponding 
eigensum contains a summation over the zeros of the combination of Bessel and Neumann 
functions. The application of the generalized Abel-Plana formula to the corresponding sum 
allows to present the Wightman function in decomposed form given by formulae (jlSp and (j23p . 
In this representations the first term on the right is the Wightman function for the wedge without 
cylindrical boundary, the term {(p{x)ip{x')) j is induced by a single shell with radius j when the 
second shell is absent, and the last terms on the right are induced by the presence of the second 
shell. For points away from the shells the last two terms are finite in the coincidence limit and 
the renormalization is needed for the first term only. By taking the coincidence limit, we have 
obtained similar representations for the VEVs of the field square and the energy-momentum 
tensor, formulae (j25p and (j34p . More symmetric decompositions are given by formulae (|3ip and 
(I42p . where the last interference term is finite everywhere including points on the shells. In the 
limit a ^ the interference parts tends to zero like a'^'^. For large values of the exterior shell 
radius, b — > oo, the interference terms in the VEVs behave as e~^'"^/6^^~^^/'^ for a massive field 
and as 1/6^"^^'^^^ for a massless one. For a wedge with small opening angle, q ^ 1, the main 
contribution into the interference parts of the VEVs comes from the summands with n = 1 and 
these parts are suppressed by the factor (a/b)'^'^. 

In section |4] we have considered the vacuum forces acting on constraining boundaries. In 
the geometry under consideration these forces are position dependent on the boundary and 
cannot be obtained by global methods using the total Casimir energy. The forces acting on 
the wedge sides are determined by the g-component of the vacuum energy-momentum tensor 
and are presented in the decomposed form (|46p . In this representation the first term on the 
right determines the force when the shells are absent and the second term is induced by the 
shells. In its turn the latter is decomposed into a single shell and second shell induced parts 
(see formula (j48p ) given by formulae (|49p . (150p . Both these forces are always attractive and do 
not depend on the curvature coupling parameter. Further we consider the forces acting on the 
cylindrical shells. These force are presented in the form ()54p where the first term on the right 
is the force acting on the cylindrical shell with radius j when the second shell is absent and 
the second term is induced by the presence of the second shell. The latter, given by formula 
([55]) . is always attractive and does not depend on the curvature coupling parameter. For large 
values of the parameter q, this part is suppressed by the factor (a/6)^'^. In the limit (po 0, 
a,b ^ oo, assuming that b — a and a(pQ are fixed, from the results of the present paper we obtain 
the corresponding formulae for the VEVs in the geometry of a rectangular waveguide. We have 
demonstrated this on the example of the interaction force between the cylindrical shells. 

Note that we have considered quantities which are well defined within the framework of 
standard renormalization procedure of quantum field theory without boundaries. We expect that 
similar results would be obtained from the model discussed in [38] where instead of externally 
imposed boundary condition the fiuctuating field is coupled to a smooth background potential 
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that reproduces the boundary condition in a Hmiting case. The generahzation of the results in 
the present paper for a scalar field with Neumann boundary conditions is straightforward. For 
this case in the expressions dU of the eigenfunctions the function cos{qn<j)) stands instead of 
sm{qn<j)) and the quantum number n takes the values 0, 1, 2, . . .. The corresponding eigenvalues 
for 7 are zeros of the function J^^ (70)1"^'^ (76) — Y^'^ (7a) 7^^(76). The formula for the summation 
over these zeros is given in [,24] . The formulae for the Wightman function and the VEV of the 
field square in Neumann case are obtained from the corresponding formulae for Dirichlet scalar 
by the replacements sm{qn(j)) cos{qn(j)), Iqn{jx) I'qnUx), Kqn{jx) K'g^{jx), j = a,b, 
and with the term n = included in the summation. In the expressions for the VEVs of the 
energy-momentum tensor this leads to the change of the sign for the second term in the figure 
braces on the right of (j35p and to the change of the sign for the off-diagonal component (136p . 
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